Three-boson problem near a narrow Feshbach resonance 
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We consider a three-boson system with resonant binary interactions and show that for sufficiently 
narrow resonances three-body observables depend only on the resonance width and the scattering 
length. The effect of narrow resonances is qualitatively different from that of wide resonances 
revealing novel physics of three-body collisions. We calculate the rate of three-body recombination to 
a weakly bound level and the atom-dimer scattering length and discuss implications for experiments 
on Bose-Einstein condensates and atom-molecule mixtures near Feshbach resonances. 
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The creation of strongly interacting atomic BEC's by 
increasing the scattering length, a, is impeded by inelas- 
tic processes |jj , the most important of which is three- 
body recombination to molecular states (see 0] , and ref- 
erences therein). For several years this process has been 
the primary motivation to study three-body physics in 
the context of ultracold gases. Recent observations of 
weakly bound dimers and their mixtures with atoms 
H,B>l3 stimulate a strong interest in atom-dimer colli- 
sions and in peculiar weakly bound trimer states (Efimov 
states) , which have been extensively studied in theory 
but have never been observed experimentally. 

In the case of a potential resonance when the scatter- 
ing length is much larger than the radius, R e , of inter- 
atomic forces, low-energy observables in a three-boson 
system are universal functions of a and the three-body 
parameter ro. In particular, for a > the rate constant 
of three-body recombination of atoms of mass to to a 
weakly bound state is given by a rec = Cha 4 /m, where 
C is a periodic function of log(o/ro) and can have any 
value between and approximately 68 0, @ ■ The atom- 
dimer scattering length is periodic in log(a/ro) and goes 
through infinity each time a new Efimov state is at the 
atom-dimer threshold 0, 0] . The three-body parameter 
ro may be determined from the three-body wavefunction 
at distances of the order of R e . For realistic interatomic 
potentials an ab initio calculation of ro is difficult and it 
is usually considered as a parameter of the theory. The 
related uncertainty makes an exact determination of the 
three-body observables impossible. Moreover, applica- 
tion of these results to three-body systems near a narrow 
Feshbach resonance should be done with care as the two- 
body scattering then depends strongly on energy. 

Feshbach resonances occur when the energy, E, of a 
pair of colliding atoms in the open channel is close to 
the energy, E res , of a quasi-discrete molecular state in 
the closed channel [llj]. In the limit when both E and 
E res are much smaller than the spacing, D ~ hr/mRl, 
between molecular levels, and the detuning E — E res is 
sufficiently small to neglect the background scattering 



[l2j . the scattering amplitude is given by [l? 
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where q = v mE/h, a = —hj/y/rnE res , and we introduce 
the length R* related to the width of the resonance by 
R* = n/y/m-y >O0. 

In this Letter we show that the resonant dependence of 
the scattering amplitude on E is of fundamental impor- 
tance for the three-body problem and develop a method 
for its solution near a narrow resonance in the case 
R* R e , \a\ ^s> R e - The three-body recombination 
rate, the positions of Efimov states, and the atom-dimer 
scattering length are shown to be functions only of the 
two-body observables, a and R* . We identify two phys- 
ically distinct regimes of the three-body scattering. In 
the regime of small detuning, \a\ 3> R* , the physics of 
a three-body collision is the same as near a wide res- 
onance. However, ro is unambiguously determined by 
the parameter R* . The regime of intermediate detuning, 
\a\ <C R*, is qualitatively different. There are no Efimov 
states for three bosons and for a > neither the atom- 
dimer scattering length nor the three-body recombina- 
tion rate shows the periodic dependence on log(a). The 
rate constant of three-body recombination to a weakly 
bound level is proportional to aJ^R* 1 / 2 and by far ex- 
ceeds the maximum value expected from the a 4 -law Q • 
Our result for aw agrees with the experimental obser- 
vation at MIT Q. 

First, we analyze the consequences of the finite reso- 
nance width for binary collisions. Although the validity 
of Eq. JTJ does not require qR* <§C 1, the amplitude (JTJ 
formally corresponds to a truncated effective range ex- 
pansion with the range R = —2R*. In the case R* ^> R e 
this range exceeds the actual radius of interatomic forces 
and the interaction acquires unusual long-range prop- 
erties. Two examples of such resonances are the 907 
G Feshbach resonance in Na (A B » 1 G, R e « 45A, 
R* w 260A) Q and the 1007.4 G resonance in 87 Rb 
(A B re 0.17 G, R e ps 85A, R* ps 320A) [HII3- Thc 
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amplitude is also a feature of optically induced reso- 
nances. A pair of free atoms can be coupled to an elec- 
tronically excited molecular level or > through a two- 
photon Raman transition reducing the spontaneous emis- 
sion, to a vibrational state of the ground-state potential 
[l8| . In contrast to magnetic resonances, the width of an 
optical one can be modified by changing laser intensities. 

The poles of the scattering amplitude (JJJ in the com- 
plex energy plane determine the positions of weakly 
bound molecular states. From the quadratic equation, 
the inverse size of the bound state is 



y/m\E b \/h = (y/l + 4R*/a- 1)/2R* 



(2) 



A true bound level exists only for a > 0. In the regime of 
small detuning we retrieve the well-known result K rs 1 /a. 
In the regime of intermediate detuning, k ~ 1/y/aR* 
and the binding energy is close to E res — —h 2 /maR*. 
Remarkably, due to the inequality R* ^> R e this level 
remains weakly bound even for a ~ R e . 

Outside the radius of the potential the normalized 
wavefunction of the bound state is given by 



(j) b (r) = (1 + 2 k R*)- 1 / 2 v / k/2tt exp(-Kr)/r, (3) 

where the preexponential factor is determined by the 
residue of the scattering amplitude at E = Ej, 13]. The 
normalization integral calculated with the wavefunction 
© tends to 1 in the regime of small detuning. We con- 
clude that atoms in the dimer are well separated and 
spend most of their time in the open channel. As the 
bound state becomes deeper the normalization integral 
decreases that reflects increasing occupation of the quasi- 
stationary level. In the regime of intermediate detuning 
the weakly bound state is this molecular level with a 
small admixture of the open-channel wavefunction. 

To solve the three-body problem we employ the zero- 
range approximation valid if the binary collision ener- 
gies are small compared to D. The key idea is to solve 
the equation for free motion but with the Bethe-Peierls 
boundary condition on the wavefunction for vanishing 
distance, r, between two atoms: —(rip)' /rip — > a~ 1 (E c ), 
where the energy-dependent scattering length defined as 



(E) 



-R*mE/h should be evaluated at the col- 



lision energy, E c , of these two atoms. This boundary con- 
dition, which can be also rewritten as ip oc 1/r— l/d(E c ), 
leads to the scattering amplitude (QJ. 

In the center of mass reference frame three bosons with 
total energy E are described by the equation 



£]*(x,y) = 



(4) 



where y is the distance between two bosons, \/3x/2 is the 
distance between their center of mass and the third atom, 
and m = h = 1. Due to the bosonic symmetry is in- 
variant under the transformations: {x — > (±\/3y — x)/2, 
y — > (V3x± y)/2}, and {x — > x, y — > — y}. Therefore, 



three boundary conditions for y — > and for \/3x±y — > 
are dependent and we need to ensure the proper behavior 
of W only for y — > 0: 



*(x,y)«[7 / - 1 -a- 1 (S c )]/(x)/4 7 r. 



(5) 



The function / contains the information about the rel- 
ative motion of the third atom with respect to the first 
two when they are on top of each other. The quantity 
E c is the total energy E excluding the energy of this 
relative motion. Accordingly, E c can be presented as a 
differential operator acting on /: 



E c = lim =E + lim = E 



f ■ 



(6) 



since a' ^(E c ) = a- 1 + R*E C , Eq. © reduces to 

* « [y- 1 - a- 1 -R*(E + V^)] /(x)/4tt. (7) 
Let us now consider the following solution of Eq. Q : 



*(x,y)=*o + J d 3 r7(r')[G £ (v^~0 



y 



Y^Ge V(x-r'/2)2 + (y±V3r'/2)2 ], (8) 



where ( x j y) is a properly symmetrized and finite 
solution of Eq. without singularities and Ge is 
the Green function of Eq. (J1J given by Ge(X) = 
—EKi(\J — EX)/8tt 3 X 2 . For negative energies K2 is an 
exponentially decaying Bessel function and for E > we 
use the convention \J — E = —iyE. A direct examination 
shows that the wavefunction (JSJ is properly symmetrized. 
In order to ensure that (jHJ reduces to Q at y — ^ we 
add and subtract from Eq. (JHJ) an auxiliary quantity 



/(x) j G £ (v/(x-r') 2 +y 2 )dV = /( x ) e -^V47ry. 

As y tends to zero we see that the terms proportional to 
1/y in Eqs. J7J| and © are equal and matching the next 
(regular) terms yields the equation for the function /: 

(-iTV 2 + L E - a- 1 + V^E - R*E)f(r) = F (r), (9) 

where Fq(t) = 47r5'o(r, 0) and Le is given by 



L E f(v) = Air 



J {G B (| r -r'|)[/(r)-/(r')] 



- 2 G E (Vr 2 + r' 2 + rr')/(r')}d 3 /. (10) 

The operator on the left hand side of Eq. © conserves 
angular momentum, and we can expand / in spherical 
harmonics to work only with a set of uncoupled equations 
for functions of a single variable r. In momentum space 
/(k) = J /(r) exp(ikr)d 3 r and Le is given by 



W(k) = (V-E + k 2 -V^E)f(k) 
2 f /(q)d 3 g 



VStt 2 J k 2 +q 2 +kq-3£:/4' 



(11) 
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Equation © was first obtained in Ref. [19| for R* = 0. 
In that case its solution is not unique |20J - at small 
distances it behaves as / ~ r _1 sinlog(r/r ) with arbi- 
trary r . According to Eq. © the corresponding colli- 
sion energy increases as \E C \ ~ 1/r 2 at small r, and for 
narrow resonances the two-body collisions would become 
off- resonant at r ~ i?*. We show below that the solu- 
tion in this case is smooth and unique. Equations 
and (JHJ determine the three-body wavefunction and all 
low-energy properties of the system. 

Consider the problem of an atom scattering by a 
weakly bound dimer (a > 0) at low collision ener- 
gies when the scattering amplitude is dominated by 
the s-wave contribution. The corresponding scattering 
length, a a d, is determined by the long-range asymp- 
tote of the atom-dimer wavefunction at zero collision 
energy. Namely, in the region y < 1/k <C x the to- 
tal three-body wavefunction factorizes into ^(xjy) w 
<f>b(y)(l — 2a a d/V3x). A comparison of this expression 
with Eqs. © and JSJ) gives f(r) oc 1 — 2a a d/V~3r at large 
r. Thus, we can find a a d by solving Eq. with E = Eb- 

(-R*V 2 T + L Eb )f(r) = 0. (12) 

To obtain Eq. ljT2*)l from Eq. JJJJ we note that — a -1 + 
yT-E h - R*E b = and * = for E < 0. 

At distances r <C 1/k we can approximate the Green 
function in Eq. JTHJ by G Eb {X) « 1/4tt 3 A 4 = G (A). 
This gives w i - The operator L has the property 
that Lor 1 ' = \{v)r v ~ l for -3 < Re (v) < 1. The function 
A(^) has two roots, v± = — 1 ± iso, where so ~ 1.00624. 
Thus, the solution of Eq. {T2J) for i?* = and r < 1/k is 

/(r) oc r" 1 sin[s ln(r/ r )] . (13) 

In the regime of small detuning Eq. 1)1 3J) remains valid 
at distances ii* « r « a, where the Laplacian in Eq. i|12|) 
can be neglected. However, for r <C i2* the differential 
operator is dominant and / should tend to a finite value 
[21j . Then, the integral operator acts as a Coulomb po- 
tential L ss A(0)/r, where A(0) = -4/V3. Note that the 
two-body collision energy \E C \ < l/R*R e <C D, which 
justifies the use of the zero-range method. 

Expression 1)1 3JI and scaling properties of Eq. 1 )12(1 
at small distances show that the long-range asymptote 
of / does not change if R* is multiplied or divided 
by exp(7r/so) « 22.7, provided that R* remains much 
smaller than a. Obviously, the atom-dimer scattering 
length has the same property. Numerical integration of 
Eq. O in the limit R* < a yields r ~0.5R*. 

In the regime of intermediate detuning the integral op- 
erator LE b in Eq. (fl"2*)l is a perturbation to the differential 
one with the small parameter l/(R*n) w ^Ja/R*. In the 
zeroth order / = const, i.e. there is no scattering. The 



first two corrections yield 

a ad /a «-8/3[l + (ll/6)v/^Rl. (14) 

In Fig. ^ we plot a a d/a as a function of R*/a calculated 
numerically from Eq. 1)12)1 . The resonances in a a d occur 
when there is a Efimov state of zero energy relative to the 
atom-molecule continuum. We find that the first Efimov 
state arises at R* /a « 2.2, the second - at R* /a w 0.07. 
With high accuracy a a d/a is periodic on a logarithmic 
scale as we further decrease R*/a. At large R*/a we 
recover the analytic result 1)14)1 . 
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FIG. 1: The ratio a a d/a versus R* /a. Dashed curve is given 
by Eq. 1)14)1 and corresponds to the limit a <C R* . Dotted 
lines indicate the points R* /a ~ 2.2 and R* /a ~ 0.07, where 
the first and the second Efimov states emerge. 

We now turn to the problem of three-body recombina- 
tion to a weakly bound level. This process occurs when 
three atoms approach each other to distances of the or- 
der of 1/k. Therefore, at energies E <C \Eb\ the recom- 
bination probability is energy independent and we can 
consider the zero-energy limit of Eq. Let the func- 
tion ^o describe three noncondensed ideal bosons in a 
unit volume. Then in the region important for the re- 
combination ^o ~ 1/V6 and Eq. 10 acquires the form 

(-i?*V2+i - a - 1 )/(r) =4tt/V6. (15) 

The solution of Eq. 1)15)1 at large r contains a term 
Aexp(iKr)/r. By virtue of Eqs. JSJ and J5J the 
corresponding part of the three-body wavefunction is 
Vf^x, y) w Ay/ (1 + 2kRF)J^kk <j>b{y) exp(inx) / x and de- 
scribes an atom and a dimer flying apart. Given the 
factor A, the derivation of the rate constant of three- 
body recombination a rcc is straightforward. We solve 
Eq. 1)15)) in momentum representation and determine A 
by the residue of the function f(k) at the pole k = K. 

In Fig. we plot a rcc /a 4 as a function of R* /a calcu- 
lated numerically from Eq. 1)15)1 . In the regime of small 
detuning our approach leads to the periodic dependence 
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of a ICC /a 4 on log(o) with the maximum value « 68 con- 
sistent with Refs. 0,0- In fact, our results for a ac i and 
for a rec in this limit coincide with those obtained by using 
the effective field theory 0, . The three-body parame- 
ter, A*, introduced there is related to R* by A* sa 6.6/i?*. 

In the regime of intermediate detuning the operator 
i?*V^ in Eq. (|f 5|l dominates and we calculate a re c per- 
turbatively. Expanding in y/ a/R*, the first two terms 
are 



a rcc m f92v / 37rVa 7 i?* 



47T 



- 1 



(16) 



Our results agree with measurements at MIT |lj . In par- 
ticular, a rcc /a 4 3> 68 in this regime. 



weakly bound level at a « 3.3R* can be used to produce 
a long-lived strongly interacting atomic gas. 
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tion for Basic Studies. 
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FIG. 2: The ratio a rcc /a 4 versus R*/a. Dashed curve is 
given by Eq. 1161 . Arrows indicate values of R*/a at which 
Qf rcc vanishes. Scatter points correspond to experimental data 
taken near the 907 G Feshbach resonance in Na at MIT [J. 

Finally, we point out that on the magnetic field axis 
the region of intermediate detuning is a factor of R* / R e 
wider than the region of small detuning and is, therefore, 
more accessible. Note, that in both MIT 1] and MPIQ 
|l5| experiments on narrow resonances a < R* . 

The method presented here can be generalized to de- 
scribe fermions and fermion-boson mixtures with differ- 
ent masses. The significant dependence of two- and three- 
body observables on R* suggests that properties of a 
many-body system should also depend on R* . 

The atom-dimer scattering can be probed by a direct 
collision of atomic and molecular clouds. It is remarkable 
that near the wide resonance associated with the first Efi- 
mov state (a 0.45i?*) the atom-dimer scattering length 
can be tuned to be an order of magnitude larger than a. 
In a degenerate atom-molecule mixture one has either a 
phase separation or a collapse depending on the sign of 
a a d- We may also consider the formation of trimers by 
adiabatic sweep across this resonance. Their detection 
can be based on the Stern-Gerlach technique 0, 0, Q 
as atoms, dimers and trimers all have different magnetic 
moments. The vanishing of the recombination to the 
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